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CN ■ Abstract 



A theoretical explanation of the observed splittings among the P states of 
charmonium is given with the use of a nonsingular potential model for heavy 



1^ . quarkonia. We also show that the recently observed mass difference between 



the center of gravity of the ^Pj states and the ^Pi state of cc does not pro- 
vide a direct test of the color hyperfine interaction in heavy quarkonia. Our 
theoretical value for the mass of the ^Pi state is in agreement with the exper- 
imental result, and its El transition width is 341.8 keV. The mass of the rj'^ 

state is predicted to be 3622.3 MeV. 
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I. INTRODUCTION 



A quantum-chromodynamic potential model was proposed by us in 1982, which 
not only yielded results for the cc and bb energy levels and their spin splittings in good 
agreement with the existing experimental data but its predictions were also confirmed by 
later experiments at the Cornell Electron Storage Ring An essential feature of our 
model was the inclusion of the one-loop radiative corrections to the quark- ant iquark 
potential, which had been derived by us in an earlier investigation Subsequently, the 
model was improved by using relativisitic kinematics [^] and a nonsingular form of the 
quarkonium potential As shown by us, in addition to the energy levels of cc and bb, our 
model also yields results in good agreement with the experimental data for the leptonic 
and El transition widths. It was further shown by Zhang, Sebastian, and Grotch |^ that 
the Ml transition widths for cc and bb obtained from our model are in better agreement 
with the experimental data than those predicted using other potential models. 

Recently the mass of the ^Pi state of charmonium has been determined by the E760 
collaboration in pp annihilations at Fermilab, and the splitting between the center of 
gravity of the ^Pj states and the ^Pi state, denoted as AMp, is found to be approximately 
—0.9 MeV. This experimental result has created much interest since it provides a new test 
for the potential models for heavy quarkonia. 

If the spin- dependent forces in the quarkonium potential could be treated 
perturbatively, the AMp splitting would arise solely from the spin-spin (color hyperfine) 
interaction. However, the spin- dependent forces are known to be quite large and, as 
observed by Lichtenberg and Potting ||^, the contributions of the spin-orbit and tensor 
interactions to AMp cannot be ignored in a nonperturbative treatment. We shall analyze 
this complex situation with the use of our model which avoids the use of an illegitimate 
perturbative treatment, and provide an explanation for the observed splittings of the 
charmonium P states. 

Several authors P-pT]] have recently shown that a theoretical value for AMp in close 
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agreement with the experimental value can be readily obtained from the spin-spin 
interaction terms in the quarkonium potential. However, since they have employed an 
illegitimate perturbative treatment, the significance of this simple interpretation remains 
an open question. 

Only a quarkonium model which is in good overall agreement with the experimental 
data can be taken seriously. Our model for heavy quarkonia satisfies this requirement. 



II. cc SPECTRUM 

Our model is based on the Hamiltonian 

H = Ho + Vp + (1) 

where 

Ho = 2{m' + pY' (2) 

is the relativistic kinetic energy term, and Vp and are nonsingular quasistatic 
perturbative and confining potentials, which are given in the Appendix. We found a trial 
wave function introduced by Jacobs, Olsson, and Suchyta |jl2[ particularly suitable for 
obtaining the quarkonium energy levels and wave functions. 

Our results for the energy-level splittings as well as the individual energy levels of cc 
are given in Tables | and |T[ For experimental data we have generally relied on the Particle 
Data Group []T3[, but for the rj^ state we have used the new result announced by the E760 



collaboration ||l4l. The two sets of theoretical results in these tables correspond to the 



scalar-exchange and the scalar-vector-exchange forms of the confining potential, given by 

K = Vs, (3a) 

and 

V,= {l-B)Vs + BVv, (3b) 
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respectively. The results obtained with the scalar-exchange confining potential are 
unsatisfactory, while the scalar-vector-exchange results are in surprisingly close agreement 
with the experimental data, including the observed mass of the ^Pi state and the AMp 
splitting. The scalar- vector mixing parameter B is found to be about ^. 

In Table fT|, we display the contributions to AMp from the various types of terms in 
the Hamiltonian (|l]) with the confining potential (|3b|) . The table shows comparable 
contributions to AMp from several sources, which brings out the complexity of this 
splitting when spin-dependent potential terms are included in the unperturbed 
Hamiltonian. The AMp splitting, therefore, does not provide a direct test of the spin-spin 
interaction in heavy quarkonia. 

In Tables IV and V, we give the results for the leptonic and El transition widths 
corresponding to the scalar-vector-exchange confining potential by using the formulae 



and 



r,,CSi^e+e-) = —-^\^{o)ni-^], (4) 
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TEiCPj^'S,) = -aelk'j\rM\ (5) 

The photon energies for the El transition widths have been obtained from the energy 
difference of the initial and the final cc states by taking into account the recoil correction. 



Our results are in good agreement with the available experimental data |T3[, and our 
prediction for the El transition width of I'Pi I'Sq is 341.8 keV. 



III. CONCLUSION 

We conclude with explanatory remarks concerning some features of our quarkonium 
potential. 



4 



A. Renormalization scheme 



We have used the Gupta- Radford (GR) renormahzation scheme [|T5| for the one-loop 
radiative corrections to the quarkonium potential rather than the modified 
minimal-subtraction (MS) scheme. The GR scheme is a simplified momentum-space 
subtraction scheme, and the parameter fi can be interpreted as representing the 
momentum scale of the physical process. This scheme also has the desirable feature that it 
satisfies the decoupling theorem [^. On the other hand, in the MS scheme /x appears as a 
mathematical parameter, and in this scheme decoupling-theorem-violating terms are 
simply ignored. 

The one-loop radiative corrections in the GR scheme can be converted into those in 



the MS scheme by means of the relation ||T5| 



a, = a. 



' ds /49 10 2 ^, m^^ 
1 + — ni + - > n — 



(6) 



47r V 3 9 3 „^ /i' 

where dg refers to the MS scheme, and ni and Uh are the numbers of light and heavy quark 
flavors. If we drop the decoupling-theorem-violating terms that appear in the MS scheme, 
we can put rii = nj and = 0, and (|) reduces to 



d, /49 10 



1 + — Uf 

47r V 3 9 



(7) 



B. Quasistatic potential 

In an earlier investigation ||^, we arrived at the surprising conclusion that while the 
quasistatic form of the quarkonium potential yields results in good agreement with the 
experimental data, this is not the case for the momentum-dependent form. This conclusion 
has also been confirmed by the recent investigations of Gara et al. and Lucha et al. 



It appears to us that the success of the quasistatic potential is related to the 
phenomenon of quark confinement. Since a rigorous treatment of quark confinement does 
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not exist at this time, we shall only offer a plausible argument. It was argued earlier |T9 
with the use of a renormalization-group-improved quantum-chromodynamic treatment 
that quark confinement can be understood as a consequence of the fact that quarks and 
antiquarks are unable to exchange low-momentum gluons. Moreover, since for the 
quark-antiquark scattering in the center-of-mass system 

where 

k = p' - p, s = p' + p, (9) 

it follows that if is allowed to take only large values, can be treated as small. This 
may be regarded as a justification for the quasistatic approximation in which terms of 
second and higher orders in s are ignored. 

Our quarkonium perturbative and confining potentials are not only quasistatic but 
also nonsingular. In the momentum space, these potentials are obtained by first expanding 
in powers of p^/(m^ + p^), and then approximating p^ as |k^. The perturbative potential 
in powers of p^/(m^ + p^) includes, among others, terms of the form 

2\ a + 6 Si ■ S2 . V 

/ P = — — — , (10) 
+ p^ 

which becomes in the quasistatic approximation 

/(k^) = ■ (11) 



+ ^k^ 



It has been observed by Grotch, Sebastian, and Zhang ||TT| that while the contribution of 
/(p^) vanishes for the P states due to the vanishing of the wave function at the origin, 
/(k^) yields a small but nonvanishing contribution for these states. Consequently, for P 
and higher angular-momentum states it would be more accurate to drop terms of the form 
(p!0|) than to convert them into the approximate form (|Tl]). We agree with the observation 
of Grotch et al. Accordingly, in the treatment of states with / 7^ we shall drop terms of 
the form ([Tl|) in the momentum-space potentials and the corresponding terms of the form 



7rr 



(12) 



in the coordinate-space potentials. 

C. Confining potential 

In our theoretical treatment, our aim has been to avoid phenomenology except in the 
choice of the long-range confining potential, which cannot be derived sufficiently accurately 
by any known theoretical technique. It is indeed remarkable that the results obtained from 
our field-theoretical perturbative potential supplemented with a phenomenological 
confining potential are in excellent over-all agreement with the experimental data including 
the AMp splitting. It should be noted that we have neglected effect of coupling of the 
energy levels to virtual decay channels and possibly other small effects. Such effects 
presumably have also been taken into account in our phenomenological confining potential. 
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APPENDIX: NONSINGULAR QUARKONIUM POTENTIALS 



The nonsingular quarkonium potentials can be obtained by appropriate 
modifications of tlie singular potentials in the momentum space, and transforming them to 
the coordinate space. The nonsingular potentials obtained by this procedure are given 
below. Some unwanted terms for states with I ^ have been dropped as explained in 
Sec. III. 



1. Perturbative quant um-chromo dynamic potential 



The perturbative potential VC, consists of the direct potential VL and the annihilation 



potential V", and in the momentum space 



W = K(k) + K"(k), 



(Al) 



where 



v;'(k) 



3k2 



^(33-2n.)ln(^^ 



3(k2 + 4m2) 



1287ra, Si -82 



-327ra 



9 k2 + 4m2 
■ (k X p) 



3«0 



k2(k2 + 4m2) 



35as 



Ois , N . / k' 



ll|(33-2,v)lng' 



+ 



3 

21a, 



In 



ISvr 



127r ) ^ \m^^ 



647ra, Si-k Sz-k - ik2Si ■ S2 



k2(k2 + 4m2 



1 + 



+ 



4a, 
37r 
3a., 



^(33-2n.)ln( 



27r 



In 



m 



(A2) 



— |^^^(l-ln2) fSi-Ss-- 
3(k2 + 4m2) ^ ^ V 4 



In the coordinate space, the potential takes the form 



(A3) 



W = KW + K"(r), 



(A4) 
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where 

v;{v) 



+ 



4as 



3q; 



3r 



27r 



-2mr , 



Gtt 



+ — (33 - 2n/) ln(^r)e-''"'' + E+{2 



mr 



la. 



+^Si • S2 ( 1 



9r 



ln(2mr)e-^"''' - K 
35a 



.(2mr)]| 



127r 



21a. 



47r 

1 - 
2a. 



|^ln(mr)e 



-2mr 



-2mr 



+ E. 



a 



+ ^(33-2n/) ln(/xr)e 

DTT 



-2mr 



E^(2mr) 



-(2mr)]| 



TT 



1 1 rv \ a 

- — ^ /i(2mr) + -^(33 - 2n^) [/i(2mr) ln(/.r) + yi(2mr)] 

ioTT / DTT 

[/i(2mr) In(mr) + gi{2mr) 



+^St Ul + ^) /2(2mr) + ^(33 - 2n/) [/2(2mr) ln(/.r) + ^2(2mr)] 



3r 



37r 



3a... 



TT 



[/2(2mr) In(mr) + 5^2 (2mr)] > , 



8a^e~^™'' / 1 



Note that the tensor operator is defined as 

5't = 3 «Ti f <T2 f — CTi • «T2, 

the functions E± are expressible in terms of the exponential-integral function Ei as 

1 



and 



/i 

/2 
^2 



e^Ei(-a;)±e ^Ei(a;) =F e ""Inx, 



1 - {l + x)e- 



(A5) 
(A6) 

(A7) 

(A8) 



1 - (1 + X + |a;2) e" 



x^ 



(A9) 



7E-K(a;)-a;^_(a;)] 



7e - 



(1 + ia;2) E+(a;) - xE_{x) 



x-^ 
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2. Phenomenological confining potential 



The scalar-exchange and the vector-exchange confining potentials in the momentum 
space are 

"1 2iS ■ (k X p; 



k4 k4(k2 + 4m2) 



(AlO) 



and 



Vv{k) = -Sir A 



1 l + |Si-S2 6iS-(kxp) Si-kSs-k-ik^Si-S 



_k4 k2(k2 + 4m2) k4(k2 + 4m2) 
The coordinate-space potentials are given by 

A 



+ 4- 



k4(k2 + 4m2) 



• (All) 



Vs r) = Ar 



2m2? 



L-S[l-2/i(2mr)], 



(A12) 



and 



-2mr 



2m2r V 3 



3A 
2'w?r 



L-S[l-2/i(2mr)] 



(A13) 



It is understood that the confining potential also contains an additive 
phenomenological constant C. 
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TABLE 1. cc energy level splittings in MeV. Theoretical splittings and parameters correspond 
to the scalar-exchange and the scalar-vector-exchange forms of the confining potential. The ex- 



perimental value of the ■0' 


— r]Q splitting 


I is not used for the determination of the cc parameters 


because of the uncertainty 


regarding the 


rj'e mass. 






Scalar 


Scalar-vector 


Expt. 


V'' - J/V' 


587.7 


588.9 


589.07=b0.13 


J/lp-Vc 


105.1 


109.0 


109 03±3 1 


i^' -v'c 


60.5 


63.5 




Xcog - J/tP 


430.5 


428.6 


428.35±1 


Xc2 - Xcl 


28.6 


44.6 


45.64±0.18 


Xcl — XcO 


80.0 


95.8 


95.43±1 


Xcog he 


-5.8 


-0.9 


-0.93±0.19±0.2 


rUc (GeV) 


2.375 


2.208 




H (GeV) 


3.329 


2.580 




as 


0.295 


0.313 




A (GeV2) 


0.183 


0.181 




B 




0.245 
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TARTF TT nn 


tjiicigy i6VGio 111 iviG V corrG&poiiQiiig, to T/Iie 


&Cd.icir-t?XCIlclIlg6 dllQ Xllc 


bCcxlcLI - VcCXOI -cXCilcLIlg 


,e loims 01 xne connnmg poxentiai. 






Scalar Scalar-vector 


Expt. 


l^Si (J/^) 


3096.9 3096.9 


3096.93±0.09 




2991.8 2987.9 


2987.9±3.1 




3684.6 3685.8 


3686.0±0.1 




3624.1 3622.3 




l'^2 (Xc2) 


3549.0 3556.0 


3556.17±0.13 


I'j'l (Xcl) 


3520.4 3511.3 


3510.53±0.12 


l=^Po(Xco) 


3440.4 3415.5 


3415.1±1 


I'PI (M 


3533.2 3526.3 


3526.2±0.15±0.2 


TABLE III. Contributions to the Xcog ~ splitting in MeV from the various types of terms 


in the cc Hamiltonian. The spin-independent, spin-spin, spin-orbit, and tensor potential terms are 


denoted as V57, V55, Vls, and Vr, respectively. 




Hamiltonian term 




Xcog - ^^"1 


^0 




19.4 


Vsi 




-9.6 


Vss 




5.2 


Vls 




-13.7 


Vt 




-2.2 



Total -0.9 
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TABLE IV. cc leptonic widths in keV. 



State 


Tee (theory) 


Tee (expt.) 




6.68 


5.36±0.29 


23^1 


3.25 


2.14±0.21 



TABLE V. El transition widths for cc in keV. The matrix elements |r/j| for these transitions 
are given in GeV~^. 

Transition J |r/j| r^i (theory) r^i (expt.) 

2=^51 ^ l^Pj 2 2.19 24.2 21.7±3.3 

1 1.96 28.1 24.2±3.6 

1.47 18.3 25.9±3.9 
l^Pj l^Si 2 1.60 293.5 270.0±33 

1 1.62 225.2 240.0±41 
1.62 105.2 92.4±42 

l^Pi ^ l^^o 1.39 341.8 



15 



